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Abstract 

Based on first-principles simulations, the structural stability and magnetic uniformity of 
Pti 3 nanoparticles encapsulated in a NaY zeolite were investigated. Among 50 stable isomers 
in the gas phase, only 15 could be accommodated into a zeolite pore and severe structural 
rearrangements occured depending on whether the solid angle at the Pt vertex bound to the 
supercage was larger than 2 sr (i.e. icosahedron). When van der Waals forces were included, 
the global minimum was found to be a new L-shaped cubic wire that is unstable in the gas 
phase. The total magnetization of the encapsulated Pti 3 decreases due to the stabilization of 
less coordinated isomers, with the majority of clusters charaterized by a total magnetization of 
2 Pb, while the majority of free clusters exhibit a threefold value. 

Within the last few decades, the number of studies related to metallic nanoparticles (NPs) 
has dramatically increased because of their applicability in catalysis, magnetic and optical de¬ 
vices and nanomedicine. 12 One of the main challenges in nanoscience is addressing the structure- 

property relationship of a structurally uniform sample of metallic NPs in a well defined size 
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range. 34 Among transition metals, encapsulated Pt-nanostructures are playing a central role in 
the fabrication of electrochemical devices, 5 including lithium batteriesj 6 because they drastically 
enhance their catalytic performance for various chemical reactions of industrial interest. 7 8 Fur¬ 
thermore, Pt-nanoystems, as clusters, nanowires, and thin films, exhibit a magnetic behavior, 912 
as a consequence of their broken translation symmetry and reduced coordination number has been 
observed. 13 Recent SQUID measurements of Pt-NPs consisting of 13±2 atoms and encaged in 
faujasite-Na zeolites show that they have a high spin state between 5 = 2 and 5 = 3, equivalent 
to a total magnetic moment of 4-6 pg, 1415 whereas XMCD studies have indicated a peak total 
magnetization of 5.9 /ig. 16 17 In any event, the sample is not uniformly magnetic; only a small 
fraction of the nominal Pt-load, specifically 15-20% as measured by XMCD spectroscopy and 
30% as measured by SQUID magnetometry are magnetic. 3D electron tomography has provided 
insights into the local structure and shape of Pt-clusters encapsulated in zeolite materials, 18 19 but 
the same authors have suggested that a local analysis, as opposed to average characterization, is 
needed to elucidate the structure-performance relationship. On this regard, atomistic simulation is 
an extremely powerful tool. Theoretical studies of metallic NPs in porous frameworks such as ze¬ 
olites, are scarce and are usually limited to one or a few atoms and to few geometrical shapes. 20 " 23 
Recent density-functional theory (DFT) calculations have suggested a Pt-icosahedron (Ih) as the 
shape grown inside zeolite pores. 24 Nevertheless, in the gas phase, this geometry is energetically 
unfavorable with respect to the buckled pyramid (BP), 25 the triangular prism (TP)J- 26 and the low- 
symmetry shape (LOW) proposed by Piotrowski et al., 27 which is the global minimum. Thus, 
a complete structural description of Ptn-NPs inserted into nanoporous materials is missing but 
highly desirable because it will provide insights into how the embedding affects the structural 
stability and then the magnetic nature of the sample. 

In this Letter, the physical properties of a variety of encaged Pti 3 isomers are calculated by 
DFT simulations, including Van der Waals dispersion forces (DFT-D). As paradigmatic example 
of porous material, a 1:4 Al/Si faujasite supercage, H3o042Al6Si24, containing two Al atoms per 
oxygen six-ring sub-structure is used as a matrix. For the first time, we demonstrate that Pt- 
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NPs@NaY are anchored to the supercage through only one atom and that their structural stability 
depends strongly on the solid angle at that anchor. When the solid angle is greater than 2 sr, i.e. the 
Ih, the isomer undergoes severe structural changes. At DFT-D, a new concave structure delimited 
by several square facets is the optimal geometry, although it is unstable in the gas phase. Among 
the ensemble of more than 15 isomers that can be inserted into a zeolite pore, the number of highly 
coordinated and highly magnetic (S>3) clusters decreases from the 70% in the gas phase to 30% 
after the embedding, due to the stabilization of low coordinated shapes. 

Model and Results 

Our numerical approach includes four different steps: accurate sampling of the potential energy 
surface of Pt-NPs in the gas phase with a focus on shapes with an high coordinate number (CN 
> 5.3) throughout an iterative metadynamics (iMT); refinement of their structural and magnetic 
properties within a DFT framework; selection and embedding within a NaY-zeolite pore; and the 
characterization of the structural and magnetic properties, including van der Waals forces (DFT-D). 

Metadynamics (MT) is an algorithm used to accelerate rare events, where the ionic dynamics is 
biased by a history-dependent potential built on a set of reaction coordinates called collective vari¬ 
ables (CVs). 28 Commonly used to explore free-energy surfaces, MT has been recently proposed 
for sampling the energy landscape of small clusters. 29 To ensure a detailed sampling of shapes in 
the gas phase, an iterative procedure (iMT) employing the coordination mumber (CN) as unique 
CV is used as in Ref. 30 The first starting configurations are well-known geometries, such as LOW, 
BP, truncated bi-pyramid (TBP), TP, icosahedron (Ih), cubooctahedron (Co) and decahedron (Dh). 
When a morphology belonging to a new structural basin is obtained, this is quenched to its nearest 
local minimum and then re-used as starting configuration for a new simulation, in an iterative man¬ 
ner. To speed up the structural search, a tight-binding second moment interatomic potetial, able to 
reproduce the geometrical details of Pt clusters, is applied. Anyway, the relative energy stability of 
each isomer is calculated according to the Broyden-Fletcher-Goldfarb-Shanno (BFGS) algorithm 
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after ionic relaxation, within a scalar-relativistic spin-polarized DFT framework, as available in 
the Quantum Espresso package. 31 Scalar-relativistic spin-polarized calculations with a non-linear 
core corrections and a small as 0.03eV Marzari-Vanderbilt smearing are able to reproduce correctly 
the geometrical features and the magnetic moment of Pti 3 with respect to relativistic approaches, 
including spin-orbit coupling. 27 32 In the case of free NPs, the simulation box contains at least 
13 A of vacuum to avoid any spurious interaction between periodic images. While for encaged 
NPs in the H3o042Al6Si24 cage, a cubic cell of approximately 28 A is used and a full BFGS ionic 
relaxation is applied to Pt atoms and the three inward O atoms of the six-ring substructure. Both 
standard DFT and corrected DFT for van der Waals forces (DFT-D), where a semi-empirical damp- 
dispersion term is introduced, 33 are performed. The NaY supercage does not exhibit any magnetic 
behavior. More details are reported in the Supporting Material pages SI-S3. 

We first discuss the results for the gas phase. An exhaustive sample of 49 different local minima 
was obtained by the EP-iMT search followed by the DFT ionic relaxation. In [TJ these isomers 
are plotted as a function of their average coordination number (CN) and maximum pair distance 
(R max)- They are labeled as "F", which represents "free objects", followed by a progressive number 
from 0 to 48. The color scheme refers to their relative stability, A E, with respect to the total energy 
of the global minimum (FO), which has been confirmed to be the low symmetry Piotrowski’s 
shape. 27 

The geometrical features of Pt-NPs embedded in NaY-zeolite, as reported in Refs.,® 16 19 can 
be summarized as follows:(i) The coordination number (CN) reaches a peak of 5.6, accordingly 
to analysis based on Fourier transform of the EXAFS ^(k) signal; (ii) the chemisorption site is on 
top of the oxygen six-ring site (SII); and (iii) the diameter should fit the dimensions of the cavity, 
and therefore be between 6 and 8 A, and have a rather spherical shape. The red box in [I] identifies 
the 14 isomers, FI to F14, that satisfy those geometrical constraints. A maximum pair distance 
of 7.5 A represents the longest distance that can be accommodated in the supercage pore, taking 
into account the Pt-0 bond length and O and Pt van der Waals radii. For the sake of completeness, 
FO will be considered in our analysis despite having a very low CN. Notably, the other low-energy 
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isomers, such as the two triangular prisms (F28 and F29), and the buckled pyramid (FI8) exhibit 
CNs that are excessively low with respect to the experimental data and thus will be not considered 
further. 



Figure 1: Energy stability of Pti 3 isomers with respect to the global minimum (FO) plotted as a 
function of their coordination number, CN and maximum pair distance R max in A. The red box 
features isomers satisfying geometrical requirements. 


All the selected clusters are reported in the first column of [2] where are grouped into their geo¬ 
metrical families. 30 The family of incomplete decahedron of 23 atoms, iDh 23 , has four isomers 
without fivefold vertex -the TBP (FI), F3, and the stellated half-star shape, F4- and only one with 
a fivefold vertex (F13). A three-fold symmetric triaugmented triangular prism (TTP), tagged F2, 
and the Dada-TTP (F6), belong to the same family. Three isomers -F5, F7, and F12- are incom¬ 
plete double icosahedra of 19 atoms, idlhi 9 . Three helicoidal shapes -F8, F10, and FI 1-, a bi-layer 
geometry (F9) -the only fully crystallographic local minimum-, and the Ihi 3 (F14) complete the 
considered ensemble of Pti 3 . Notably, only FI has a A E below 0.5 eV while all the other isomers, 
but Ih, lying between 0.6 and 1.5 eV. With respect to magnetic properties, the selected free Pto 
clusters show a magnetic character. Their total magnetization (TM), calculated as the difference 
between the majority and minority charge densities integrated over the real space, is reported in 
|2} The F0-F14 sample has a TM average of 5.6 /is, with roughly the 45% with a TM of 6 /i«. 
Notably, wire-like minima such as F41 and F42, which arise from the compenetration of octahedra 
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of 6 atoms, have a very high TM (10-12 [1b). Unfortunately, these shapes are too long (R max 
A,) to be grown into the supercage. 
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Figure 2: Selected Pti 3 isomers, classified into their geometrical families. The left columns rep¬ 
resents the gas phase (F prefixes), the middle column represents Pti 3 inserted into the zeolite su¬ 
percage (Z prefixes); and the third column represents refinement at the DFT-D level (D prefixes). 
For each cluster, the total magnetization in Bohr magneton (/is), the relative energy with respect 
the global minimum (in eV), and the solid angle at the anchor are listed below their identification 
label. The coordination of the anchor itself is listed in brackets. 
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The F0-F14 isomers were then inserted into the NaY zeolite at the SII supercage site, and ionically 
relaxed. We first note that the chemisorption is always through only one Pt vertex anchored to only 
one O atom of the six-ring of the supercage, in vert good agreement with experimental data. As 
depicted in|2j ZO to Z14 identify the embedded clusters in zeolite using a standard DFT framework, 
whereas DO to D14 refer to the case where London dispersion corrections are included. A label 
R indicates isomers that reconstruct significantly from their starting configuration. The relative 
stability of each minimum is calculated with respect to the global minimum, ZO-c and DOR-c, at 
DFT and DFT-D, respectively. DOR-c, which is an L-shaped motif delimited by six squares, arises 
from ZO-c when dispersion forces are introduced. In [3} free and encaged minima are plotted as 
functions of their CN and TM, and the color scheme refers to their energy stability with respect 
to the global minimum. At the DFT level, the overall coordination of Pt isomers, and thus their 
magnetic properties, is almost preserved after the isomers are embedded with the exception of the 
unstable geometries that undergo a major reconstruction, i.e., Z6R, Z7R, Z8R, Z13R, and Z14R, 
as in [2 
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Figure 3: Average coordination number (CN) with respect to the total magnetization, TM (j Ub). 
The color refers to energy stability of Ptl3 -NPs in the gas phase (left), embedded in an NaY 
zeolite (middle) and in zeolite with dispersion forces DFT-D (right). A E (eV), relative energy 
stability, to the global minima F0, ZO-c and DOR-c, respectively. The region between the red lines 
identifies the experimental range of coordination. 
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Within a DFT-D framework, the net effect of adding dispersion forces is a decrease in the nearest- 
neighbor average distance, a preference for assuming a cage-like shape associated with a higher 
binding energy, and the shift towards lower coordinated isomers is made evident in [3} Compared 
with the DFT case, the distance between the anchor and the center of the oxygen ring is always 
shorter by 0.1-0.3 A . Simultaneously, all the embedded NPs tend to move toward the center of the 
cage, where the average distance between the cluster centers of mass and the geometrical center of 
the pore is 0.2-0.3 A . The dispersion forces thus induce a few local atomic rearrangements that are 
sometimes significant. The most relevant case is the global minimum ZO-c. It evolves toward an L- 
shaped structure with six square facets, DOR-c, linked to the supercage through a three-fold vertex. 
Notably, DOR-c is unstable in the gas phase and folds again into the original LOW shape. However, 
its counterpart chemisorbed with a three-fold vertex, DO-b, remains geometrically stable, although 
it loses 0.4 eV in energy. At the DFT-D level, the A E is spread up to 2.5 eV and it is strongly 
correlated to the coordination number where A E > 1.5 eV are for CN > 4.7. 


Discussion and Conclusion 

To fully understand the geometrical stability of zeolite-embedded Pt-NPs, we investigate the role 
of the vertex anchoring the cluster to the supercage. At the DFT level, the chemisorption through 
a 3- and 4-fold coordinated vertex usually leads to negligible changes in the overall cluster shape, 
whereas chemisorption through a 5-fold coordinated vertex leads to major structural transforma¬ 
tions. To be more quantitative, the encapsulation energy E enc and the solid angle at the anchor, Cl, 
were calculated. E enc is the energy gain due to the embedding, where values greater than 1.3 eV 
indicate structural rearrangement of the whole cluster. The values of Cl reported per each isomer 
in [2] represents the angles subtended by the polygon obtained by connecting all the Pt atoms coor¬ 
dinated to the vertex, as projected onto the unit sphere around it. Our simulations show that only 
solid angles Cl smaller than 0.8 sr are fully stable, while values greater than 2 sr are completely 
unstable from a geometrical point of view; i.e., the polygon underneath the vertex/anchor can be a 
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triangle or a quadrilateral, but not a rectangle or planar pentagon. Among clusters with Cl < 0.8 sr, 
the best examples are the truncated bi-pyramid (Zl) and the symmetric TTP-based structure (Z2); 
they have a three-fold rotational symmetry around the Pt-anchor, and a E enc ~ 0.9 eV. A non-planar 
pentagon, with a borderline value of Cl, is allowed, as in the case of the bilayer (Z9) which has a 
E enc ~ 0.9 eV. When Cl > 2 sr and E enc >1.4 eV, severe rearrangements occur to decrease Cl to 
approximately 1 sr, as depicted in [2j Generally speaking, structural reconstruction lowers the TM 
with the sole exception of the isomer arising from the Ih reconstruction. Notably, the Ihi 3 indeed 
evolves towards a iDh 23 upon a three-fold anchor, a distorted half-star, Z4. Not surprisingly, the 
reconstructed Ih exhibits a TM of 4 [1b similarly to Z4. A full description of those transformations 
comes as Supporting Information. Those structures that reconstructed at the DFT level because of 
the embedding continue to change toward less coordinated shapes at the DFT-D level. Notably, the 
incomplete decahedra are geometrically stable while helicoidal shapes tend to rearrange. Details 
of these transformations are provided at pages S5-S7 of the Supporting Information. 

The structural rearrangements due to the embedding and the inclusion of van der Waals forces 
strongly affect the magnetic properties of Pti 3 -NPs@NaY. Although TM ranges between 2 and 8 
Hb as in the gas phase, the average TM is about 3.9 /i# against 4.8 /is at the DFT level, and 5.6 
Hb in the gas phase, as shown in |4j The decrease in the TM is associated to the stabilization of 
low coordinated shapes. D7 is the only case which increases its TM to 8 jis due to a different 
orientation within the supercage. Thus far, at the DFT-D level, only six isomers have a TM above 
6 Hb ■ TBP (Dl); TTP (D2); ilhi 9 (D12) and the helicoidal D10-R, which have a spin state of 3; 
and the other two incomplete icosahedra D5 and D7 which show a high magnetization of 8 /l#. 
Anyway, only four of them are highly coordinated (CN>5). 


Bartolome and co-workers 16 suggested that the magnetic behavior of Pt-NPs could be influenced 
by a charge-transfer between the O atoms of the cage and the cluster. Here, we observe a charge 
transfer of approximately 0.2 e“/atom from the O toward the Pt vertex bound to the six-ring sub¬ 
structure. The involved oxygen exhibits an atomic polarizability of ± 0.4 jlB that contributes 
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Figure 4: Proportion of considered isomers in terms of their total magentization TM (Hb), in the 
gas phase (cyan), in zeolite (blue) and DFT-D zeolite (magenta). 

to the total magnetic moment of the whole system. The presence of two A1 atoms close to the 
vertex/anchor of the NP results in a Pt-O distance peaked at 2.25 A in excellent agreement with 
experiments. 14 

In conclusion, we investigated the energy landscape of Pti 3 using a metadynamics approach 
and found an exhaustive set of highly coordinated isomers with a maximum pair distance shorter 
than 7.5 A. Fifteen different isomers were selected and then embedded in a NaY zeolite supercage 
model with different orientations when possible. Each Pti 3 -NP adsorbs at the SII site with a vertex 
anchored to the inward O atom of the six oxygen ring substructure. The effect of the embedding 
was the destabilization of NPs linked to the supercage throughout solid angles larger than 2 sr, as 
occurs in the case of the icosahedron. Van der Waals forces were included to determine the stability 
of each isomer, and their magnetic properties were subsequently estimated. In the case of encap¬ 
sulated Pt-NPs, a new concave structure was found to be the lowest energy structure. As in the 
gas phase, the embedded Pt-NPs exhibit a magnetic behavior, with a total magnetization between 
2 and 8 fis- However, due to the stabilization of low coordinated isomers, the majority exhibits 
a TM of 2 Hb against a value of 6 /i# calculated in free clusters. The high magnetic moment of 
Pti 3 -NPs in zeolite, was mainly due to four highly coordinated isomers: a truncated bi-pyramid, a 
triaugmented triangular prism, and two incomplete double icosahedra. Nonetheless, these isomers 
are quite unfavorable from an energetic point of view, which explains why only a fraction of the 
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total Pt load is magnetically active in very good agreement with available experimental data. 
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